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A lattice perspective of kaon phenomenology * 
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Abstract. I review recent lattice computations of the matrix element relevant for K" — mixing and 
discuss the advantages of fermions with an exact chiral symmetry to compute K ~* tttt amplitudes. 
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1 Introduction 



Kaon decays are sources of key informations to test the 
flavour structure of the Standard Model and to search for 
new physics. Recently the experimental results for indi- 
rect CP violation in if ^ tttt decays and those from B de- 
cays and oscillations were compared to test the Cabibbo- 
Kobayashi-Maskawa picture of flavour mixing and CP vi- 
olation milOl- In these analyses the experimental result 
for e is connected to the fundamental parameters of the 
underlying electroweak theory through the lattice deter- 
mination of the relevant hadronic matrix element of the 
AS = 2 effective Hamiltonian. Its bag parameter Bk has 
been computed in the quenched approximation over the 
past decade with a remarkable accuracy. A summary esti- 
mate of Bk is given in the first part of this talk, followed 
by a review of recent computations. 

Quantitative tests of the same quality have not been 
possible so far in the AS = 1 sector for the so-called 
AI =1/2 rule in K —^ tttt decays or for the direct CP vi- 
olation parameter e'/e. In this sector a more complicated 
blend of ultraviolet and infrared effects renders the deter- 
mination of the relevant matrix elements very challenging 
and prevented so far reliable computations with standard 
regularizations. In the ultraviolet, power-divergent sub- 
tractions can be needed to construct the renormalized op- 
erators that enter the effective Hamiltonian In 
the infrared, the continuation of the theory to Euclidean 
space-time and the use of finite volumes in numerical sim- 
ulations generate a non-simple relation between the physi- 
cal amplitudes and those computed on the lattice p^lF^HTT)] . 
The use of lattice fermion discretizations which preserve 
an exact chiral symmetry at finite lattice spacing greatly 
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simplifies some of these computations and allows one to 
attack them as shown in the second part of this talk^. 

Before entering into details, it is important to stress 
that lattice QCD allows one to perform non-perturbative 
computations of masses and matrix elements from first 
principles, with control, at least in principle, over all sources 
of systematic errors. These uncertainties can be systemat- 
ically reduced by exploiting the properties of the underly- 
ing quantum field theory and/or by using more powerful 
computers, but without the introduction of any model- 
dependent assumption or parameter. Examples are given 
by the 0(a)-improved actions and operators and by the 
non-perturbative techniques used to renormalize bare ma- 
trix elements. Up to now many computations are still per- 
formed in the so-called quenched approximation, i.e. drop- 
ping the vacuum polarization effects in the Monte Carlo 
simulations. For some of them this represents the only 
source of systematic uncertainty not under control. 

2 TTTT decays 

Non-leptonic K tttt amplitudes can be parametrized as 

Is 



T[K+ - 
T[K° - 



TT TT 



(1) 



iS2 



where Si and Aj are the tttt phase shifts and the isospin 
amplitudes for / = 0, 2. Direct and indirect CP violations 



^ For lack of space only a selection of topics are reviewed 
here, chosen following my personal taste and competence. 
Please see Ref. |11| for an up-to-date report on rare kaon de- 
cays and Ref. |12| for a broader overview on the lattice activity 
in kaon physics. 
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are parametrized by 

, _ J_ j, ReA2 /ImA2 _ Im^p N 
^ ~ ReAo \ReA2 ReAo) ^ ' 

and 

^ T[Kl ^ (w)o] 

r[X5 ^ (7r7r)o] ^''^ 

respectively. Experimental results reveal = tt /2 + 62 — 
60 w 7r/4, a Z\/ = 1/2 selection rule ^ 22.2, the 

presence of direct and indirect CP violation in nature 

Re(£7e) = (16.6 ± 1.6) x 10"^ [WTIj 

|£| = (2.282 ± 0.017) X 10^3 . (4) 

The large value of |yloM2| and that of e'/e can be ex- 
plained within the Standard Model only if the strong in- 
teractions crucially affect these non-leptonic weak transi- 
tions (see 16,17,18 for recent reviews). 

2.1 K°-K° mixing: e 

By using the OPE, the AS = 2 effective Hamiltonian is 
given by 

^oi^=' = ^|f^Cr(M)ai(M)+h.c., (5) 

where the expression of the Wilson coefficient Ci(/i) as a 
function of the fundamental parameters of the underlying 
theory is known at NLO 19 and the corresponding bare 
four-ferniion operator is given by 

Oi = is-t^P.d){s-t^P.d) , P± = . (6) 

Its matrix element between and K'^ 

(K°\^^\K')^^F|,M],BK (7) 

encodes the long-distance QCD contributions to £. So far 
the most precise and careful computation of ]3k (/i) in the 
quenched approximation has been carried out by JLQCD 
PI7| (see also |2J|). They used staggered fermions, one- 
loop perturbation theory to renormalize the operator and 
degenerate down and strange quark masses. Their results 
at fixed lattice spacings are reported in Fig. ^ together 
with their continuum extrapolated value 

i3^(2GeV) = 0.63 ± 0.04 (8) 
= 0.86 ± 0.06 , (9) 

where Bk is the renormalization group-invariant bag pa- 
rameter. The final error of 0.04 is much larger than the 
statistical ones at fixed lattice spacing. The amplification 
is due to the fact that with these fermions 0{a^) effects in 
the matrix element are found to be large, and a consistent 
continuum limit of the results can be obtained only if also 
O(a^) corrections are allowed in the extrapolation fit. By 
using chiral perturbation theory, Sharpe estimates a 5% 



systematics due to the degeneracy of the quark masses 
|22j . At present the systematic error due to the quench- 
ing approximation is not under control. A crude estimate 
based on quenched chiral perturbation theory and par- 
tially quenched simulations suggests to include a further 
15% error in Eq. ^ 23j. The error on Bk can be put un- 
der control and significantly reduced only with full QCD 
simulations with realistic light-quark masses. Despite a 
great effort in the last few years in the lattice commu- 
nity j23p24j ■ the simulation of dynamical quarks is still 
very challenging and a breakthrough in numerical algo- 
rithms may still be needed to reach interesting light-quark 
masses. 

In the last year several collaborations have been comput- 
ing Bfc in quenched QCD by using sophisticated regu- 
larizations, non-perturbative renormalization techniques 
and allowing for ^ P^26II27| . The aim is to check 
the staggered result and to pin down the systematic error 
within the quenched approximation. Even if a significant 
comparison can be performed only when all systematics 
(within the quenching approximation) will be under con- 
trol also in these computations, it is interesting to analyse 
and compare the results to appreciate their potentiality. 
Two collaborations '25''26' computed Bk with Ginsparg- 
Wilson fermions 28,29,30,^ for the first time. These 
regularizations simultaneously preserve chiral and flavour 
symmetries at finite lattice spacing |32| . As a consequence 
Oi renormalizes multiplicatively and its matrix elements 
are 0(a)-improved. The plain Neuberger action [301 has 
been used in Ref. PJI for a lattice of linear extension 
L ~ 1.5 fm, a spacing a ~ 0.09 fm and for degener- 
ate light-quark masses. The RI/MOM non-perturbative 
renormalization procedure has been implemented to com- 
pute the logarithmic divergent renormalization constant. 
The result obtained is reported in Fig. ^ Even with a 
large error, it is in very good agreement with the contin- 
uum extrapolated staggered result, suggesting that also 
for this matrix element discretization effects may be mild 
with Neuberger's fermions. A NNC-HYP overlap fermion 
action has been implemented in Ref. 26 for lattices of 
linear extension L ~ 1.3 fm, a ~ 0.08, 0.12 fm and for de- 
generate light-quark masses. The results arc lower but still 
compatible within two sigmas with the continuum extrap- 
olated staggered one and are in the same ballpark as those 
obtained with domain-wall fermions with a finite fifth di- 
mension 33 , 3^ . Even if both studies require more work 
to properly assess the magnitude of the various systematic 
errors, they confirmed that Ginsparg- Wilson fermions can 
be very effective in computing phenomenological quanti- 
ties despite their numerical cost. 

In the last few years the ALPHA collaboration has been 
computing Bk with the goal of reaching a very precise 
determination of Bk within the quenched approximation. 
They are using twisted-mass fermions, which break ffavour 
symmetry and violate parity but allow for a multiplica- 
tively renormalization of Oi. A non-perturbative deter- 
mination of the renormalization constant and its running 
has been completed, while the computation of the matrix 
element is still under way. A preliminary result has been 
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Fig. 1. B^S(2GeV) versus tlie lattice spacing in units of the 
Sommer scale ro |35l from the reference computation, black 
circles I2UI (see also also black squares |21|'l. and from recent 
determinations: red diamonds 125) . green circles |26) . red tri- 
angle |27|. orange diamonds "33", orange squares )34|. 



reported in ^7] and is shown in Fig. ^ Also in this case 
the result obtained at a ~ 0.09 is in very good agreement 
with the continuum extrapolated staggered result. 

2.2 The AI ^1/2 rule 

By using the operator product expansion (OPE), the CP- 
conserving ZiiS* = 1 effective Hamiltonian above the charm 
threshold is given by P^37U38| 



ttAS=1 



Gp 

71 



where the Wilson coefficients C±(^) are known at the 
NLO prnilin] and the bare operators are 



Oh 



(10) 



c) 



The contributions that arise when the top quark is inte- 
grated out are heavily suppressed by CKM factors and 
can be neglected. 0± belong to different chiral multiplets 
and are CPS-even. In correlation functions at non-zero 
physical distance, 0± cannot mix between themselves or 
with other four-fermion operators if the regularization pre- 
serves chiral symmetry, but only with the dimension-six 
operator [51C] 



(ml 



md{sP-f-d) + ms{sP^d) 



The renormalized operators are 



(11) 



(12) 



where Z± (/i) are logarithmic-divergent renormalization con- 
stants and bj!: are suppressed by a factor as- No power- 
divergent subtractions are needed to renormalize 0± when 
fermions with an exact chiral symmetry are used 0. 
For rus ^ m^, 



Qrn = {ml - ml)df. 



Asd 



rud 



■Vf. 



rrid 



(13) 



rrid 



and it does not contribute to matrix elements that pre- 
serve four- momentum 4,5 . 

If the charm is integrated out, not only potentially 
large contributions of 0{p? /ml) are neglected, but ultra- 
violet power divergences can arise in the renormalization 
pattern of the relevant four-fermion operators. In this case 
the AS = 1 effective Hamiltonian can be written as 



rTAS=l 
-"eff 



10 

1=1 



The so-called QCD-penguin operators are 

S3,5 = (S7p^-rf) E (^^M^T?) 

q—u,d,s 

q—u,d.s 



(14) 

(15) 
(16) 



(see Refs. |16II17| for definitions of the other operators). 
At non-zero physical distance, mixing with two lower- 
dimensional operators 

Qp = md{sP+d) + ms{sP-d) (17) 
Qa = rnd{sFi^,^a^j_yP+d) + ms{sF^,yai^,yP^d) 

can occur and power-divergent subtractions are needed 
even with Ginsparg- Wilson fermions. 

With Wilson fermions, only CPS and flavour symme- 
try can be used to determine the renormalization pattern 
of 0±. Even with an active charm, a quadratic diver- 
gent contribution needs to be subtracted in the parity- 
conserving sector 



PC 



3 



where 



0£^(a^) = 2±(m) o 



+ blQr 



Qs = {mu - mc)sd 

Qt = {rriu - mc)sai^,„F^^d 



(18) 



(19) 
(20) 



and Oj_ are four-fermion operators with wrong chirality. 
In this case only the flavour part of the GIM mechanism 
survives thanks to the explicit breaking of chiral symmetry 
It is interesting to notice that power divergences can 
be mitigated with twisted-mass Wilson fermions |41| . 
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At leading order in chiral perturbation theory, the ma- 
trix elements needed for the AI ~ 1/2 rule can be ex- 
tracted from three-point correlation functions, thus avoid- 
ing the infrared problem that affects the direct computa- 
tion of the K TTTT matrix elements on the lattice 0]. 
No large cancellations among leading order terms are ex- 
pected in the ratio IA0/A2I; an enhancement should there- 
fore be visible already at this order. As a result, a com- 
bined use of fermions with an exact chiral symmetry and 
chiral perturbation theory can be the starting point to 
attack the AI =1/2 rule. 

In the past years the RBC and the CP-PACS collab- 
orations have studied the AI —1/2 rule^ with domain- 
wall fermions with a finite fifth dimension |33II42| . They 
have computed the if — > tt and K ^ matrix elements 
for the operators of the AS — 1 effective Hamiltonian in 
Eq. H14|l and used LO chiral perturbation theory to recover 
the physical amplitudes. Albeit with large statistical and 
systematic errors, both groups demonstrated that a con- 
trolled numerical signal can be obtained for these matrix 
elements. The systematics of these important results can 
be put under control by using fermions with an exact chi- 
ral symmetry, lighter quark masses, and by considering the 
effective Hamiltonian with a dynamical charm. A numer- 
ical feasibility study with Neuberger's fermions is under 
way 031 ■ It is also conceivable that the relevant low-energy 
constants of the weak chiral Lagrangian can be extracted 
by studying the weak interactions in the e-regime [4311441 

For direct CP violation, both the ultraviolet and the 
infrared problems are more severe. An active charm does 
not mitigate the ultraviolet renormalization, and diver- 
gent power subtractions are necessary. The cancellation 
between two large competing contributions from Qq and 
Qs renders s' /s very sensitive to higher order corrections 
in chiral perturbation theory. Leading-order terms may 
not be sufficient to reach a reliable prediction in the Stan- 
dard Model gZlllHl. 

I warmly thank A. Vladikas and T. DeGrand for send- 
ing their results and for interesting discussions about their 
work. Many thanks to L. Lellouch for a critical reading of 
the manuscript and to the organizers for the very stim- 
ulating atmosphere of the conference. The participation 
to the conference was supported in part by the Euro- 
pean Community's Human Potential Programme under 
contract HPRN-CT-2000-00145, Hadrons/Lattice QCD. 
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